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Abstract 

This study investigates the effectiveness of ventilation design strategies for general hospital wards in terms of 
virus removal capacity. A typical semi-enclosed six-bed general ward of Hong Kong hospitals and three 
respiratory viruses, namely Middle East respiratory syndrome coronavirus (MERS-CoV), severe acute respira- 
tory syndrome coronavirus (SARS-CoV) and H1N1 influenza virus, were chosen for the computational fluid 
dynamics (CFD) simulation of airflow field and virus dispersion inside the ward. The results demonstrated that 
the location of an infected patient would affect the infection risks to other occupants and healthcare workers 
inside the same hospital ward, and an increased air change rate in the ward could reduce the risk of infection 
from direct contact and inhalation. It was found that an air change rate of 9h‘ could effectively minimize the 
deposition and floating time of respiratory virus particles while maximizing energy efficiency. This study 
should provide a useful source of reference for the hospital management to mitigate the risk of infection 


with MERS or other airborne transmitted viruses through better ventilation design strategies. 
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Introduction 


According to a statistical report by the Hong Kong 
Hospital Authority,' the in-patient discharges and 
deaths were continuously increasing from 2003 to 
2013. The 2012/2013 overall number of in-patient dis- 
charges and deaths was 1,027,005, and that of day- 
patient was 516,127 in Hong Kong. Among all patients, 
33.4% were day-patients. 

To prevent nosocomial or healthcare-associated 
infections, especially airborne ones, hospital hygiene 
and infection control are necessary. The Centres for 
Disease Control and Prevention (CDC) provides guid- 
ance to help healthcare personnel to follow standard, 
contact, and airborne precautions when caring for hos- 
pitalized patients with known or suspected viral infec- 
tions.* Effective prevention measures, e.g. an airborne 
infection isolation room (AIIR), are especially crucial 
for control of acute respiratory infectious threats. 


Proper ventilation also plays a key role in infection 
control by minimizing airborne bacteria and viruses. 
Both the spread of severe acute respiratory syndrome 
coronavirus (SARS-CoV) during the largest nosoco- 
mial SARS outbreak in Hong Kong and the recent out- 
break of Middle East respiratory syndrome (MERS) in 
the South Korean hospitals revealed that airborne dis- 
ease transmission through inefficient hospital ward ven- 
tilation systems can lead to dire health consequences.* > 
For a balanced ventilation that delivers indoor air 
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quality and energy efficiency, general hospital wards 
have been designed to meet certain air change require- 
ments.”*° '° For instance, an air change rate (ach) ran- 
ging from 2 to 6h’ is suggested to help decrease local 
mean age of air,'! while 4h~' is recommended for 
energy savings.'*'? However, as new information 
becomes available, current air change requirements 
for hospital wards should be re-evaluated and updated. 

This study investigates the effectiveness of ventila- 
tion design strategies for general hospital wards in 
terms of virus removal capacity. A computational 
fluid dynamics (CFD) simulation of a typical general 
ward of Hong Kong hospitals was conducted and three 
respiratory viruses, namely MERS-CoV, SARS-CoV 
and HIN1 influenza virus, were chosen. The findings 
can be used by the hospital management to minimize 
cross-infection risk while reducing energy consumption 
of ventilation. 


Bioaerosol drag force 


The motions of spherical and non-spherical bioaerosols 
in a ventilated space can be calculated using equation 
(1) by integrating the force balance on the bioaerosols 
in terms of the drag force per unit particle mass per 
relative velocity Fp (N s kg~' m~'), where g is gravita- 
tional acceleration (m s °), F, is the additional acceler- 
ation force per unit particle mass (N kg~') if any, vj 
and v, are the velocities of the virus and air (m s_'), 
respectively, j, is the molecular viscosity of air 
(kgm~'s~'), d, is the equivalent bioaerosol diameter 
(um), Re, is the Reynolds number for bioaerosols in 
an airflow field, , is the air density (kg m ) and pp is 
the virus density (=1100 kg m~+).'* 


dv, —~ Pa 
b= P=) 4 802 Pa) Gos 
dt Pb 1) 
F 18a se CpRe, Re Padp|Vo = Val ( 
p=—=_ *—5, 2) ha = 
a} pp 24 Ha 


The drag coefficient Cp for the bioaerosols is defined 
by equation (2) 


Cp = —; Re, < 1 (2) 


Equations (1) and (2) are used in the CFD simula- 
tion to determine the motions of droplet nuclei under a 
Lagrangian scheme. 

The drag constant for the bioaerosols Kp in equation 
(2) is given by equation (3) 


RSS 


Kp = (3) 
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Validity of equation (3) was established for equivalent 
bioaerosol diameters d, = 0.69—-6.9 zm over a bioaerosol 
size range and further examined for Bacteriophage Phi 
X174 (ATCC 13706-B1) with d; = 0.054 um. A combined 
experimental and simulation method described by Wong 
et al.'* was performed. Before the experiment, agar plates 
were prepared in such a way that the entire preparation 
was covered by the host bacterium Escherichia coli 
(ATCC 13706). Plaque assays for bacteriophages were 
then conducted. The number of bacteriophages collected 
was measured by plate counting.'” Figure 1 shows that 
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Figure 1. Bioaerosol drag constant Kp against equivalent 
bioaerosol diameter d,. 
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Figure 2. Bioaerosol drag constants and absolute errors for 
Bacteriophage Phi X174. 
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equation (3) can be extrapolated for predicting the disper- 
sion of smaller sized bioaerosols (down to d, = 0.054 wm). 
When using Kp=0.001458 (Figure 2), the deviation 
between simulated and measured fractional counts 
inside a ventilated chamber was insignificant (p > 0.99, 
paired ¢-test). The corresponding fractional bias (FB) «, 
and normalized mean square error (NMSE) «> for 
Kp= 0.001458 were 9.2% and —3%, respectively, as illu- 
strated in Figure 2.'° 


Using equation (4), where the projected image area, 
length and width are A,, /; and /, respectively, the equiva- 
lent bioaerosol diameter d, can be determined from elec- 
tron micrographs.'”'* Electron micrographs of samples 
of MERS-CoV, SARS-CoV, HINI1 influenza virus and 
Bacteriophage Phi X174 are shown in Figure 3.'? *” 


|A max (I, 1 
dp =2 Vaspect = 2) (4) 


min (1, 1) 


CFD simulation 


A typical semi-enclosed six-bed general ward cubicle 
(7.5m (L)x6m (W)x 2.7m (A)) employed for the 
CFD simulation is shown in Figure 4. As illustrated 
in Figure 5, there were four ceiling air inlets with 
mounted diffusers in the cubicle, and the corridor was 
the air-outlet. Four different air change rates (ie. 
ach = 3, 6, 9 and 13h7') were used, while the tempera- 
ture and relative humidity (RH) of the supply air were 
285K and 80-95%, respectively.7* Six probable viral 
emission points, representing the positions of six 
patients lying on their respective beds (i.e. Man 1, 
Man 2,..., Man 6), were chosen as the source locations. 
The viruses investigated were MERS-CoV, SARS-CoV 
and HIN1 influenza virus. The metabolic rate of a 
reclining patient was assumed to be 0.8 MET (ce. 
46.6W m~).** Half of the heat (23.3W m~) trans- 
ferred from patient skin surface by convection was 
assumed.*° Patient beds were set as rectangular box 


Figure 3. Electron micrograph of reference viruses. (a) MERS-CoV'? (b) SARS-CoV”? (c) HINI influenza virus”! and (d) 


Bacteriophage Phi X174.7° 
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CORRIDOR 


WINDOWS 


Figure 4. A six-bed general ward cubicle (dimensions in mm). 


for the worst scenarios that airflow under the beds was 
fully blocked by medical equipment and relation instal- 
lations. This arrangement may cause higher numbers of 
bioaerosol particle deposited on other patients. 

The simulation of airflow field and virus dispersion 
inside the ward cubicle was done through the use of 
CFD software FLUENT 14. An Eulerian—Lagrangian 
framework was used to solve the gas—solid two-phase 
flow problem, i.e. an Eulerian scheme for the prediction 
of a steady-state airflow field of the ward cubicle, fol- 
lowed by a Lagrangian approach for the determination 
of virus particle movements. 

In the Eulerian framework, a continuous phase of 
the induced airflow field was obtained with a second- 
order solution scheme. The renormalization group 


(RNG) k-e turbulence model was tested to be an appro- 
priate choice among the RANS turbulence models, and 
therefore, it was adopted to determine the air turbu- 
lence in the field. The model offered better accuracy 
and stability in cases of low Reynolds number and 
near-wall flows, and this model was found suitable for 
indoor airflow simulations also.*°?’ To couple with the 
pressure and velocity fields, the pressure implicit with 
splitting of operator (PISO) algorithm was employed, 
and the convection term was discretized using a second- 
order upwind scheme. 

To optimize the simulation quality and speed, three 
reference grid sizes namely fine (i.e. 2618 k), moderate 
(i.e. 1143k) and coarse (i.e. 1071 k) were constructed 
for determining a suitable mesh size. The fine gird 
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Figure 5. CFD configurations of a six-bed general ward cubicle. 


size was with a mesh skewness less than 0.25 to ensure 
an excellent cell quality. The moderate and coarse grid 
sizes were doubled and quadrupled by a second-order 
method for increasing the simulation speed. To assess 
the mesh quality using linear grid stretching, the 
asymptotic range of convergence Cgsym, applied 
was based on the grid convergence indexes (GCIs) 
with a relative error of computed average mass flow 
rate Epms.* 

If the value of Cgsymp in equation (5), where f; is the 
safety factor, r, is the refinement ratio and c, is the 
theoretical order of convergence, is approximately 1, 
then the grid quality is fine. 


GCTfine = evel > GCTeoarse — ane : 
. ne re — 
(5) 
G Cl eoarse 


Casymp = TAL 
(GCI fine)1; 


The GCI analysis results from the simulation were 
51%, 37% and 1.187 for GCIeoarses GCTjine ANA Casymps 
respectively. 

In the CFD simulation, the transmission pathway 
of virus-laden respiratory droplets expelled by 
sneezing was predicted (i.e. d;, initial virus-laden 
respiratory droplets diameter = 8.3 um, v,=50m s |, 
pp=1100kg m=? and n,, amount of virus particles 
expelled by sneezing=10,000 virus particles).'*°? +! 
General patient wards are recommended at an RH 


30-60%.° Within a short period of time (<0.1s) after 
emission, the droplets would evaporate to droplet 
nuclei.** These nuclei were the dried-out residual of 
droplets possibly containing infectious pathogens.* 
Only a small proportion (<10%) of the virus-laden 
droplets of a total number of 73,000—1,000,000 droplets 
expelled by a vigorous sneeze was assumed; neither 
aggregation with other particles nor cluster of virus 
would be formed from the dried-out nuclei at such 
low concentration.***° Lipid-enveloped human corona- 
virus 299 E would remain alive by a half-life of 67h at 
an RH of 50% and an air temperature of 20°C.*° The 
survival rate of the three virus was expected to be 100% 
in simulation times less than 100s. Other details for 
simulations and boundary conditions are summarized 
in Table 1. 

Table 2 exhibits the equivalent bioaerosol diameters 
d, for the virus droplet nuclei examined in this study. 
A one-way coupling was applied in the prediction to 
prevent the effect of particles on the continuum airflow. 
Each virus particle was tracked separately for its pos- 
ition and velocity by a discrete phase model (DPM). 
A previous study confirmed that the isotropic dis- 
crete random walk (DRW) model was effective and 
accurate in modelling bioaerosol dispersion and distri- 
bution due to turbulent fluctuations in the flow.*’ For a 
coagulation effect of bioaerosol particles in this study, a 
very low volume fraction (<3000cm~*) was kept in the 
ward cubicle to reduce collisions of the virus particles in 
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a turbulent flow.***? Both the amounts of virus par- 
ticles exhausted to the corridor n, and deposited on 
the ward room surfaces ny were counted.*° A perfect 
sink boundary condition was applied to the ward room 
surfaces in order that the virus particles impinged onto 


Table 1. CFD simulations and boundary conditions. 


7.5m (L) x 6m (W) x 2.7m (A), 
RNG k-e turbulence model, 
standard wall function. 

1,143,766 cells 

1.24kg s~! (for ach =3), 2.48kg 
s-! (for ach=6), 3.72kg s~! (for 
ach=9) and 5.37kg s_! (for 
ach = 13), 285K (air 
temperature) 

0.31 kg s~' (for ach =3), 0.62kg 
s! (for ach= 6), 0.93 kg s~! (for 
ach=9) and 1.34kg s_! (for 
ach = 13), 285K (air 
temperature) 


Computational domain 
Mesh configuration 


Total supply airflow rate 


Each inlet airflow rate 
(0.6m x 0.6m) 


Four diffuser 
(0.6m x 0.6m) 


Four-way spread-type, supply jets 
had an angle of 15° from ceiling, 
adiabatic and reflect boundary 
type. 

Pressure-outlet, 295 K (backflow 
temperature), adiabatic, escape 
boundary type. 

No slip wall boundary, 295K 


Corridor (6m x 2.7m) 


Walls, ceiling, floor and 


beds (surface temperature), adia- 
batic, trap boundary type. 
Patients No slip wall boundary, 23.3 W 


m ~*, trap boundary. 


Patient mouths 
(0.05 m x 0.05 m) 


Single-shot release with an upward 
exhalation velocity of v,=50m 
s_', initial virus-laden respira- 
tory droplet diameter 
d;=8.3 um, n,= 10,000 virus 
particles, density of bioaerosol 
particle p,=1100kg m~* 


RNG: renormalization group. 


Table 2. Virus information. 
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the solid surfaces would be perfectly trapped with no 
reflection and desorption.*! 


Deposited ratio 


Direct contact with the viruses deposited on ward room 
surfaces and inhalation of the viruses suspended in the 
air are two potential transmission routes of virus-laden 
airborne particles expelled by sneezing.** The deposited 
ratio ry 1s selected as a measure for evaluating the con- 
tact transmission; the movement of sneeze particles can 
be represented by equation (6), where r, is the exhausted 
ratio for particles exhausted to the corridor and r, is the 
elapsed ratio for particles suspended in the air 
r+r,try=l (6) 

Ratios r,, r, and rg are given by equation (7), where 
Ns, Ney Ng and ng are the amounts of virus particles 
expelled by sneezing, exhausted to the corridor, elapsed 
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Figure 6. Description of the bioaerosol removal process in 
a general ward. 


Equivalent Aspect Evaporation Evaporation Evaporation 

bioaerosol ratio Drag time at 0% time at 50% time at 90% 
Species ATCC diameter d,, (tum)* (aspect) constant Kp RH (s) RH (s) RH (s) 
MERS-CoV - 0.167 £0.012 1.27 0.013945 3.48 x 107° 9.55 x 107° 1.82 x 107! 
SARS-CoV - 0.1375 + 0.009 1.052 0.009453 3.48 x 10° 9.55 x 10°? 1.82 x 107! 
HINI influenza virus - 0.124 + 0.0001 1.203 0.0199 3.48 x 10-7 9.55 x 10 1.82 x 107! 
Bacteriophage Phi X174 13706-B1 0.054+ 0.014 1.012 0.001458 - - - 


MERS-CoV: Middle East respiratory syndrome coronavirus; SARS-CoV: severe acute respiratory syndrome coronavirus. 


“Standard errors shown. 
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(a) Temperature (K) (b) Air velocity vq (ms) 
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(c) Air velocity v, (m s"') 
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Figure 7. Simulation results of the ward with ach=6h7'. (a) Temperature distribution, (b) air velocity distribution and (c) 
flow pathlines from air supply inlets. 
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Figure 8. MERS-CoV pathways for six source locations with ach=6h7'. (a) Man 1, (b) Man 2, (c) Man 3, (d) Man 4, 
(e) Man 5 and (f) Man 6. 
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(with respect to air suspension) and deposited on ward 
room surfaces, respectively 


(7) 


Figure 6 graphs the ratios r,, r, and ry against time. 
The elapsed time rt, (s) required for all virus particles to be 
exhausted or to deposit on ward room surfaces can be 
determined from the figure as represented by equation (8). 


rq(Ta) = 0 (8) 
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Results and discussion 
Airflow field 


Figure 7 shows the simulated temperature, air velocity 
distribution and the flow pathlines inside the ward with 
ach=6h~'. An increasing temperature gradient was 
found from the ceiling to the floor inside the ward. As 
shown in Figure 7(a), warmer air observed near the 
patients (due to body temperature) was diffusing out to 
the corridor. The average room temperature was about 
296K. Stagnant air, with velocity less than 0.05m s~', 
was found near the window side (i.e. Wall 2) and the 
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Figure 9. Deposited and exhausted ratios of MERS-CoV for six source locations. (a) Man 1, (b) Man 2, (c) Man 3, (d) Man 


4, (e) Man 5 and (f) Man 6. 


Downloaded from ibe.sagepub.com at DREXEL UNIV LIBRARIES on June 5, 2016 


Yu et al. 


lower floor level (i.e. below 0.3m) as presented in 
Figure 7(b). The flow pathlines in Figure 7(c) confirm 
that there was insufficient ventilation to remove any 
virus particles in those areas. Generally, the airflow dir- 
ection was from the inner part of the ward to the 
corridor. 


Cross-infection from surface deposition 


Figure 8 indicates the potential risks of cross- 
infection with MERS-CoV through air pathways from 
different infected patients inside the ward. Patients stay- 
ing on the same side of an infected patient, especially 
the one located next to the corridor (i.e. Man 1 or 
Man 2), would have a higher chance of cross- 
infection. Two different virus pathway flows in the 
simulation due to the asymmetric diffuser locations 
are highlighted. Cases exhibited in Figure 8(c) and (e) 
show the virus moved along floor surface of the 
ward but in cases shown in Figure 8(d) and (f), virus 
would pass over nearby patients’ heads, then flew to 
the corridor. 

Figure 9 exhibits the values of deposited ratio ry and 
exhausted ratio r. of MERS-CoV on all surfaces inside 
the ward with ach=3, 6, 9 and 13h !. The results sug- 
gested that the deposition of virus particles was depend- 
ent on the location of an infected patient, and the virus 
particles would deposit mainly back on the source 
patient. If the infected patient was located near the cor- 
ridor (i.e. Man 1 or Man 2), the virus particles would 
likely be exhausted to the corridor. On the contrary, if 
the infected patient was located in the inner part of the 
ward (i.e. Man 5 or Man 6), the possibility that the 
virus particles would deposit on the wall surfaces or 
other patients was higher. 

Moreover, more virus particles would remain on the 
sneezing patient at lower air change rates, especially for 
ach=6h~! (rg>0.7 for all source locations), whereas 
more virus particles would deposit on the wall surfaces 
or be exhausted to the corridor at higher air change 
rates (i.e. ach=9h~' or 13h7'). Although the possibil- 
ity of cross-infection among patients and healthcare 
workers (HCWs) was lower at a higher ach (up to 
9h~'), HCWs should clean all patients and ward sur- 
faces regularly, regardless of ach. 


Cross-infection from inhalation 


Figure 10 graphs the MERS-CoV removal processes in 
the ward for the six source locations (i.e. Man 1, Man 
2,..., Man 6) and ach=3, 6, 9 and 13h7!. The results 
showed that virus particles from an infected patient 
located near the corridor were likely to be exhausted 
to the corridor. Meanwhile, lower r, (i.e. higher ry) 
and longer elapsed time t, were associated with the 


MERS-CoV particles from the sneezing of an infected 
patient located in the inner part of the ward, and that 
indicated a higher risk of MERS-CoV infection 
through direct inhalation of particles or indirect inhal- 
ation of re-suspended particles. Besides, higher r, and 
shorter t, were found with increasing air change rates. 
If ach was increased from 3 to 13h~', t, could be shor- 
tened by more than 30s, and the risk of cross-infection 
from inhalation could be effectively lowered. Similar 
dispersion and deposition results were observed for 
SARS-CoV and HIN] influenza virus particles and pic- 
torialized in Figures 11 and 12. The simulation per- 
formances demonstrated that virus particles with a 
relatively small equivalent bioaerosol diameter 
(d, ~ 0.1 wm) had similar particle dispersion and depos- 
ition characteristics in a general hospital ward. 

Figure 13 plots the average elapsed time t, against 
ach for MERS-CoV, SARS-CoV and HIN1 influenza 
virus. The results showed that t, could be significantly 
shortened by increasing the air change rate in the ward 
(R* > 0.9). However, the threshold or optimal t, in the 
ward could not be determined as the infectious doses 
(ID50) of the three viruses varied in wide ranges (i.e. 
180, 1800 and 180 virus particles for MERS-CoV, 
SARS-CoV and H1N1 influenza virus, respectively, to 
cause a 50% infection).** *° Using the ASHRAE stand- 
ard for an AIIR (i.e. ach= 12h!) as a safety measure,° 
the corresponding elapsed time was about 34s. When 
applying the ASHRAE (ach=4h~') and CIBSE 
(ach=6h') standards for a general patient room, the 
average values of t, were found to be 70s and 61s, 
respectively.°* The elapsed time doubled when ach 
dropped from 12 to 4h7', and thus doubling the poten- 
tial inhalation risk. Based on the median value in 
accordance with both ASHRAE and CIBSE standards, 
the maximum ach in a general hospital ward should be 
9h! (t,=48 s) for the needs of maximizing energy effi- 
ciency and minimizing infection risk. Furthermore, it 
should be noted for other means to minimize the pos- 
sibility of cross-contamination in hospital wards, such 
as the installation of ultraviolet germicidal irradiation 
(UVGI) lamps for the destruction of viral nucleic 
acids.*° 


Conclusion 


This study should provide a useful source of reference 
for the hospital management to mitigate the risk of 
infection with MERS or other airborne transmitted 
viruses through better ventilation design strategies. 
The results of this study demonstrated that the location 
of an infected patient would affect the infection risks to 
other occupants and HCWs inside the same hospital 
ward, and an increased air change rate in the ward 
could reduce the risk of infection from direct contact 
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Figure 10. MERS-CoV removal processes for different source locations and air change rates. 


Downloaded from ibe.sagepub.com at DREXEL UNIV LIBRARIES on June 5, 2016 


Yu et al. 


Man 1 


Ss S Ss Ss 
o iv = a oe ms 


o 
=) 
—) 
+ 
Ss 
s 
~ 
— 
s 
o 


20 40 60 80 100 


Man 2 
S Ss = J 
i—) nN > a o — 
© eee 
So N > a ~ - 


o 


20 40 60 80 100 


= 


0 2 40 60 80 100 


0 2 40 60 80 100 0 2 40 60 80 100 0 2 40 60 80 100 0 2 40 60 80 100 
1 1 1 
08 08 08 
(oa) 0.6 0.6 0.6 
i= 
S 04 04 04 
0.2 02 02 
eee 0 0 0 
0 2 40 60 80 100 0 2 40 60 80 100 0 2 40 60 80 100 0 2 40 60 80 100 
1 1 1 1 
08 08 0.8 08 
tr 06 0.6 0.6 0.6 
< 
S 04 04 04 04 
02 02 02 02 
0 as 0 0 0 


° 
8 
s 
3 
g 
S 
° 
8 
s 
Es 


80 


3 


04 


0.2 


Man 5 
cs e¢ee 
So wR > a oo -_ 
os © 
i—) a te - 


04 


0.2 


o 


20 40 «60 #80) «100 


_ 


0 2 40 60 80 100 


0 20 40 60 80 100 


11 


0 2 40 60 80 100 0 20 40 60 80 100 0 20 40 60 80 100 
1 1 1 1 
0.8 08 08 08 
OO 06 0.6 0.6 0.6 
i=} 
$ 04 04 04 04 
02 02 02 02 
0 - 0 0 0 
0 2 40 60 80 100 0 2 40 60 80 100 0 20 40 60 80 100 0 2 40 60 80 100 
ach = 3 bt! ach = 6b! ach=9 h'! ach = 13 hr! 


x-axis: Time (s) 
y-axis: Deposited ratio ra, exhausted ratio re 


Figure 11. SARS-CoV removal processes for different source locations and air change rates. 
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Figure 12. H1N1 influenza virus removal processes for different source locations and air change rates. 
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MERS: Middle East respiratory syndrome; SARS: severe acute 
respiratory syndrome; AIIR: airborne infection isolation room. 


and inhalation. For a typical semi-enclosed six-bed gen- 
eral ward of Hong Kong hospitals, an air change rate 
of 9h! could effectively minimize the deposition and 
floating time of respiratory virus particles while max- 
imizing energy efficiency. In order to minimize the pos- 
sibility of cross-contamination in hospital wards, 
installation of UVGI lamps is also recommended. 
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